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ABSTRACT 

A  common  technique,  referred  to  as  channeled  imaging  polarimetry  (CIP),  enables  the  snapshot  acquisition  of  the  2- 
dimensional  Stokes  parameters  of  an  arbitrary  scene  or  sample.  It  achieves  this  by  amplitude  modulating  the  Stokes 
parameters  onto  various  interference-based  spatial  carrier  frequencies.  While  this  technique  has  utility,  it  often  suffers 
from  low  signal-to-noise  ratios,  in  remote  sensing  scenarios,  due  to  its  requirement  of  narrow  spectral  bandwidth 
illumination  (<  3  nm  in  the  visible).  This  paper  discusses  one  hardware  implementation  that  can  be  utilized  to  overcome 
this  limitation.  Consequently,  an  overview  of  the  theoretical  and  experimental  development  of  a  unique  Sagnac 
interferometer  is  discussed.  Both  laboratory  and  outdoor  data  are  included  to  demonstrate  the  instrument's  ability  to 
measure  polarization  in  arbitrary  scenes.  Inclusion  of  blazed  diffraction  gratings  inside  the  interferometer  enables  white- 
light  interference  fringes  to  be  generated.  By  incorporating  these  gratings,  the  operational  bandwidth  of  the  interference 
fringes  can  exceed  approximately  300  nm;  two  orders  of  magnitude  greater  than  previous  CIP  implementations.  Lastly, 
by  modifying  the  diffraction  grating,  the  sensor  becomes  capable  of  snapshot  multispectral  imaging.  This  will  be  briefly 
discussed,  with  both  a  theoretical  description  and  experimental  data. 
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1.0  INTRODUCTION 

Instantaneous  acquisition  of  the  Stokes  polarization  parameters  is  of  great  interest  in  many  areas  of  remote  sensing  [  1  ].  A 
Stokes  imaging  polarimeter  is  capable  of  obtaining  either  the  partial  or  complete  polarization  state  of  a  scene  via  four 
Stokes  parameters.  These  parameters  express  the  state  of  polarization  in  a  4x1  matrix,  defined  as 
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where  v,  v  arc  spatial  coordinates  in  the  scene,  .So  is  the  total  energy  of  the  beam,  St  denotes  preference  for  linear  0“  over 
90°,  S2  for  linear  45°  over  135°,  and  Sj  for  circular  right  over  left  polarization  states.  Each  is  defined  by  an  addition  or 
subtraction  of  intensity  measurements  that  represent  different  analyzer  states,  and  complete  characterization  requires  at 
least  four  such  states  be  measured.  These  measurements  can  be  acquired  over  time  by  use  of  a  rotating  element  (division 
of  time)  with  a  single  imaging  lens  and  focal  plane  array  (FPA).  Such  an  implementation  yields  a  compact  and  relatively 
inexpensive  instrument.  However,  these  sensors  are  susceptible  to  misregistration  errors  caused  by  platform  or  scene 
motion,  since  in  many  applications.  Stokes  parameters  are  acquired  from  moving  platforms. 

One  method  of  acquiring  all  the  Stokes  parameters  simultaneously,  discovered  by  K.  Oka,  is  by  amplitude  modulating 
the  Stokes  parameters  onto  carrier  frequencies  [2].  Oka  established  that  the  complete  Stokes  vector  can  be  encoded  onto 
various  interference  fringes  with  the  use  of  two  Savart  plates  located  within  a  collimated  space  of  an  optical  system  [3]. 
This  provides  the  advantage  of  being  snapshot,  while  also  offering  inherent  image  registration  since  the  Stokes 
parameters  are  encoded  on  coincident  fringe  fields.  However,  one  previously  unresolved  concern  is  that  the  carrier 
frequency’s  visibility  decreases  as  the  coherence  length  of  the  incident  illumination  is  decreased.  Therefore,  high- 
visibility  fringes  are  generated  only  when  the  light  is  quasi-monochromatic,  leading  to  a  reduced  signal-to-noise  ratio  in 
remote  sensing  applications.  This  is  especially  a  concern  for  operation  of  the  sensor  in  the  thermal  infrared  (3-12  //m)  [4). 
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In  this  paper,  we  discuss  the  theoretical  and  experimental  development  of  a  dispersion  compensated  Sagnac 
interferometer.  In  section  2,  we  discuss  the  use  of  a  Sagnac  interferometer  as  a  replacement  for  a  Savart  plate  in  K. 
Oka’s  original  implementation.  In  section  3,  the  dispersion  compensated  Sagnac  polarization  interferometer  (DCPSI)  is 
introduced,  and  its  theoretical  background  is  provided.  In  section  4.  an  interesting  aspect  of  the  DCPSI  is  investigated,  in 
that  it  can  be  applied  for  use  as  a  multispectral  imager.  Lastly,  section  5  highlights  experimental  results  from  both 
systems. 


2.0  SAGNAC  INTERFEROMETER  AS  A  SAVART  PLATE 

A  Sagnac  interferometer  with  polarization  optics  can  duplicate  the  shearing  and  optical  properties  of  a  Savart  plate  (SP) 
[5,  6,  7],  Such  a  polarization  Sagnac  interferometer  (PSI)  is  portrayed  in  Fig.  1 . 


Fig.  1.  A  polarization  Sagnac  interferometer.  The  distance  between  the  WGBS  and  mirrors,  M,  and  VL.  is  denoted 
d\  and  d2,  respectively.  A  shear  (Spy)  is  produced  when  d\  i-  d:.  The  case  d\  >  d->  is  illustrated,  with  d,  =  d2  +  a.  [8] 

The  PSI  contains  two  mirrors,  M)  and  M2,  at  the  output  faces  of  a  wire  grid  beam  splitter  (WGBS).  It  also  contains  a  lens 
with  a  focal  length. /,bj,  which  is  used  to  image  the  source  or  scene  onto  a  focal  plane  array  (FPA).  The  two  beams  are 
separated  by  a  shear,  SPSJ,  when  dt  i-  d2.  If  dt  =  d2  +  «,  where  a  is  some  distance  that  mirror  M|  was  translated,  then  the 
shear  is 

S„,  =  -hex  (2) 

Note  that  this  shear  is  independent  of  the  wavelength  of  the  incident  illumination  [9],  The  optical  path  difference  (OPD) 
between  the  two  sheared  beams  emerging  from  the  PSI  is  given  by 

OPD  =  Srsi  sin (0)  =  Srsiff.  (3) 


When  the  two  sheared  beams  are  imaged  onto  the  FPA  by  the  lens,  they  combine  and  produce  interference  at  the  image 
plane.  The  intensity  at  the  FPA  is  described  by 
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where  S0,  Si,  and  S3  are  implicitly  dependent  upon  spatial  coordinates  x,  and  y,.  Fourier  transformation  of  this 
interference  pattern  produces  multiple  carrier  frequencies  in  the  spatial  frequency  domain.  Filtration  of  these  peaks, 
followed  by  an  inverse  Fourier  transform,  enables  the  incident  Stokes  parameters  to  be  extracted.  [3]  A  carrier  frequency. 
Upsi,  can  be  described  as 
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Notable  is  the  inverse  relationship  to  the  wavelength.  This  dispersion  in  the  spatial  carrier  frequency  causes  a  visibility 
reduction  in  the  fringes  that  arc  modulating  the  incident  Stokes  parameters.  When  the  carrier  frequency  visibility  is 
reduced  to  zero,  the  incident  Stokes  parameter  becomes  irrecoverable.  This  typically  occurs  for  spectral  bandwidths 
exceeding  approximately  3  nm.  To  remedy  concerns  related  to  the  narrow  spectral  bandwidth,  the  shear  must  be  made 
linearly  proportional  to  the  wavelength,  such  that  Sps,  «  Ay,  where  y  is  some  optical  thickness,  analogous  to  a. 


3.0  DISPERSION  COMPENSATION  IN  THE  SAGNAC  INTERFEROMETER 


Removing  the  dispersion  in  the  PSI’s  carrier  frequency  can  be  achieved  by  introducing  two  blazed  diffraction  gratings 
into  each  arm  of  the  interferometer.  Diffraction  gratings  have  been  used  in  past  experiments  to  generate  white-light 
interference  fringes,  specifically  for  optical  testing  [10].  The  new  interferometer,  referred  to  as  the  dispersion- 
compensated  polarization  Sagnac  interferometer  (DCPS1),  is  depicted  in  Fig.  2.  In  this  case,  the  WGBS  to  mirror  M|  and 
Mi  separations  (d,  and  di  in  Fig.  1 )  are  set  equal,  and  the  shear  is  generated  exclusively  by  the  diffraction  gratings,  G| 
and  Gi.  A  ray’s  diffraction  angle,  after  transmission  through  cither  grating,  is  calculated  by 

sin(0)=-^-,  (6) 

d 

where  f>  is  the  diffraction  angle  of  the  ray  as  measured  from  the  grating's  normal,  m  is  the  order  of  diffraction,  and  d  is 
the  period  of  the  grating.  Since  d  is  typically  large  (>=  30  /an),  small  angle  approximations  can  be  used  to  simplify  Eq.  6. 
yielding  9  ~  mA/d  . 


Fig.  2.  DCPSI  with  blazed  diffraction  gratings.  G,  and  G2,  positioned  at  each  output  of  the  WGBS.  Inclusion  of 
the  gratings  generates  a  shear  that  is  directly  proportional  to  the  wavelength. 

In  the  DCPSI,  the  spectrally  broadband  beam  reflected  by  the  WGBS  is  initially  diffracted  by  grating  G:  into  the  +1 
order.  These  dispersed  rays  propagate  through  the  system  to  grating  G,.  where  they  arc  diffracted  into  the  0-order.  thus 
removing  the  diffraction  angle  imposed  by  Gj.  The  rays  exit  the  system  parallel  to  the  optical  axis,  but  offset  by  a 
distance  proportional  to  >.v„,  where  x„  is  some  constant  related  to  the  DCPSI’s  configuration.  Conversely,  the  beam 
transmitted  by  the  WGBS  is  initially  incident  and  diffracted  by  Gt.  and  exits  the  system  offset  by  -Ax„.  The  functional 
form  of  the  shear,  calculated  by  unfolding  the  optical  system,  is  expressed  as 
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where  a.  b ,  and  c  represent  the  distances  between  G|  and  Mi,  Mj  and  Mi,  and  M ; and  G2,  respectively.  Consequently,  the 
DCPSI  can  generated  a  shear  that  is  linearly  dependent  upon  wavelength  to  first-order.  The  intensity  on  the  FPA  can 
then  be  expressed  as 
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where  the  total  intensity  pattern  is  a  summation  from  the  minimum  to  maximum  order  of  a  blazed  grating,  such  that  the 
maximum  achievable  order  is  (cM|)sin(;t/2).  where  /,  is  the  minimum  wavelength  passed  by  the  optical  system  and 
Su  (m),  and  S\  (m)  are  the  Stokes  parameters  weighted  by  the  diffraction  efficiency  (DE)  of  both  gratings  before 

integration  over  wavelength. 
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wavelengths  passed  by  the  optical  system.  The  carrier  frequency, 
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Therefore,  inserting  the  blazed  diffraction  gratings  into  the  optical  system  enables  cancellation  of  the  dispersion  tenn  in 
the  denominator.  Another  interesting  aspect  of  the  carrier  frequency  is  that  it  depends  upon  the  order  of  diffraction,  m, 
which  can  be  exploited  for  multispectral  imaging. 

4.0  MULTISPECTRAL  DCPSI 

In  the  polarimetric  DCPSI,  it  is  advantageous  to  utilize  a  single-order  blazed  grating  optimized  for  the  +1  diffraction 
order,  which  should  be  blazed  for  a  wavelength  in  the  middle  of  the  spectral  operating  regime.  However,  if  a  multiple- 
order  diffraction  grating  is  used  (i.e.  a  grating  with  a  larger  phase  depth),  then  larger  diffraction  orders  (eg.,  m  =  4,  5,  6 
etc.)  would  be  observed  in  the  system.  Each  one  of  these  orders  corresponds  to  a  different  blaze  wavelength,  such  that 
=  Xt/m.  Since  each  order  that  is  observed  on  the  image  plane  contains  a  different  spatial  carrier  frequency,  then  each 
diffraction  order  can  be  treated  as  a  unique  spectral  passband. 


An  illustration  of  the  ideal  diffraction  efficiencies  for  the  gratings  used  in  the  experiments  is  depicted  in  Fig.  3.  Fig.  3  (a) 
portrays  the  diffraction  efficiencies  of  the  single-order  blazed  grating  used  in  the  DCPSI,  while  Fig.  3  (b)  provides  the 
diffraction  efficiencies  for  the  multiple-order  blazed  grating  used  in  the  MSI. 


ideal  muoretical  diffract,  or  efficiency,  singie-ordec  grating  «ea:  theoretical  diffract  ion  efficiency.  mJliple-oroer  grating 


Fig.  3.  (a)  Diffraction  efficiency  (percent)  for  the  blazed  gratings  used  in  the  experimental  results  of  the  (a)  single- 
order  DCPSI  and  (b)  multiple-order  MSI.  [II] 

The  configuration  for  the  MSI  is  depicted  in  Fig.  4.  and  is  similar  to  the  DCPSI  per  Fig.  2  except  for  the  inclusion  of  the 
multiple-order  gratings.  Additionally,  a  linear  polarizer  (LP)  at  45°  is  also  included  to  uniformly  polarize  the  incident 
illumination.  This  makes  the  intensity  pattern  on  the  FPA  equal  to. 
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where 


$(*)  =  J  D£2(/l,/w)(S0(,l)-i-S,(/l))rf/l  (14) 

and  the  dominant  orders  experimentally  observed  in  the  system  correspond  to  the  ceiling  (Ce)  of /|//mm.  where  A|  is  the 
first  order  blaze  wavelength  of  the  diffraction  grating.  This  changes  the  maximum  limit  of  the  summation  from  dl),m in  to 
Ce[z i observed  in  Eq.  13.  For  small  incident  S2,  the  measurement  corresponds  to  S0,  or  the  unpolarized  component 
of  the  incident  illumination.  It  should  be  noted  that  a  polarization  insensitive  design  can  be  produced  by  replacing  the 
WGBS  with  a  standard  50/50  beamsplitter. 


Fig.  4.  MSI  configuration.  The  single-order  gratings  of  the  DCPSI  are  replaced  with  multiple-order  (higher  phase- 
depth)  gratings.  [II] 


5.0  EXPERIMENTAL  RESULTS 


To  verify  the  operating  principles  of  the  DCPS1  and  MSI,  an  experimental  setup  was  configured  for  each  variant  of  the 
system.  For  the  DCPS1.  the  polarimetric  accuracy  was  assessed,  while  with  the  MSI,  the  spectral  reflectance  accuracy 
was  studied. 

5.1  DCPS1 

The  experimental  configuration  for  the  DCPSI  is  depicted  in  Fig.  5.  From  Fig.  5,  spatially  and  temporally  incoherent 
light  is  setup  through  the  use  of  a  reflective  diffuser.  A  dichroic  polymer  linear  polarizer,  oriented  at  45°,  followed  by  an 
achromatic  quarter-wave  retarder  (QWR)  at  #G,  is  used  to  generate  input  polarization  states  for  testing  and  validation. 
The  WGBS  has  a  clear  aperture  of  approximately  21  mm  and  is  anti-reflection  coated  for  400-700  nm.  Mirrors  M|  and 
M;  consist  of  1/1 0^  wave  optical  flats,  while  a  collimating  lens  (fc  =  50  mm)  is  used  to  image  the  object  plane  to  infinity. 
An  objective  lens  (Job/  =  200  mm)  is  used  for  re-imaging  onto  the  FPA,  thus  yielding  a  field  of  view  (FOV)  of 
approximately  +/-  1°.  The  IR  blocking  filter  rejects  wavelengths  beyond  750  nm  from  entering  the  objective  lens. 


Fig.  5.  Experimental  setup  for  laboratory  testing  of  the  DCPSI  using  white  light.  The  distances  a,  b.  and  c  arc  14.4 
mm,  57.2  mm,  and  12.9  mm,  respectively. 

First,  an  experimental  verification  of  the  white-light  nature  of  the  carrier  frequency  was  preformed.  This  consisted  of 
placing  a  monochrometer  in  front  of  the  instrument.  Scanning  the  monochrometer  from  460-700  nm  in  10  nm 
increments,  while  measuring  the  frequency  of  the  interference  fringe,  yielded  the  results  depicted  in  Fig.  6  for  the 
DCPSI’s  carrier  frequency  ( Uncpsi )•  For  comparison,  the  calculated  carrier  frequency,  as  a  function  of  wavelength,  for  a 
standard  PS1  (£/«■,)  is  also  included.  Using  the  theoretical  data,  the  carrier  frequency's  total  peak-to-peak  variation  for 
the  uncompensated  PS1  was  calculated  to  be  1 1.31E3  m  '.  Conversely,  the  total  variation  in  the  carrier  frequency  of  the 
DCPSI  was  measured  to  be  23.61  m1.  Consequently,  the  total  variation  in  the  carrier  frequency  for  the  DCPSI  is  over 
two  orders  of  magnitude  lower  than  in  the  uncompensated  PSI. 
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Fig.  6.  Measured  carrier  frequency  vs.  wavelength  for  the  DCPSI  (solid  dark-gray  line)  and  the  theoretical  carrier 
frequency  for  a  PS1  (dotted  light-gray  line).  The  PSl's  carrier  frequency  was  set  to  equal  the  DCPSI’s  carrier 
frequency  at  a  wavelength  of  600  nm. 

Data  of  outdoor  targets  were  also  obtained  with  the  instrument.  Data  taken  of  a  moving  vehicle  is  depicted  in  Fig.  7. 
Here,  the  achromatic  quarter-wave  retarder  (QWR)  was  used  at  the  entrance  of  the  system,  positioned  in  front  of  the 
WGBS,  enabling  St  in  Eq.  8  to  be  converted  into  S2.  Thus,  a  complete  linear  polarimeter  that  can  simultaneously 
measure  So,  Si,  and  S2  is  produced  [12]. 


Fig.  7.  (a)  Raw  data  of  a  vehicle  outdoors,  (right  panel)  Processed  polarization  data  of  the  vehicle,  calculated  from 
the  data  in  (a).  From  the  image  in  the  first  row  and  first  column  and  moving  clockwise:  Sn,  DOLP,  S2fS(h  and  5|/S<i- 

5.2  MSI 

The  experimental  configuration  for  the  MSI  is  depicted  in  Fig.  8.  Here,  spatially  and  temporally  incoherent  illumination 
is  used  to  illuminate  a  target  that  is  to  be  measured  in  reflectance.  Furthermore,  the  objective  lens’  focal  length  is 
increased  to  =  250  mm,  such  that  the  higher  diffraction  orders  appear  as  a  lower  frequency.  Otherwise,  all  other 
specifications  for  the  WGBS  and  mirrors  arc  identical  to  the  DCPSl’s  experimental  configuration. 


Fig.  8.  Experimental  setup  of  the  MSI  with  a,  h,  and  c  set  to  13.3  mm,  51.4  mm,  and  19.1  mm,  respectively.  An  IR 
blocking  filter  keeps  wavelengths  >  750  nm  from  entering  the  system  [II] 

In  order  to  measure  the  relative  spectral  responsivity  of  the  MSI's  spectral  passbands,  the  exit  slit  of  a  radiometrically 
calibrated  monochrometer  was  placed  within  the  object  plane.  The  output  of  the  monochrometer  is  depicted  in  Fig.  9  (a). 
This  enabled  the  carrier  frequencies  to  be  observed  and  measured,  where  the  absolute  value  of  their  amplitude  in  the 
Fourier  domain  was  utilized  to  determine  their  spectral  response.  These  responsivity  measurements  arc  depicted  in  Fig.  9 
(b).  As  predicted  from  the  ideal  diffraction  efficiencies  depicted  previous  in  Fig.  3  (b).  four  bands  arc  present  within  the 
system. 


Monocfrromotef  Output  t*  Wnvotength  Retativo  Spectral  Response  vs  Wavelength 


Fig.  9.  (a)  Measured  monochromctcr  output  in  W/m2.  (b)  Measured  relative  spectral  response  of  the  different 
passbands  within  the  MSI.  Also  included  (solid  light-gray  line)  is  the  relative  response  of  the  FPA  multiplied  by 
the  measured  transmission  of  the  IR  blocking  filter.  [  1 1] 

In  order  to  investigate  the  accuracy  of  the  multispectral  data,  a  reflectance  measurement  from  a  healthy  and  unhealthy 
leaf  was  taken  with  the  MSI.  and  compared  to  a  calibrated  non-imaging  spectrometer.  The  measured  passbands, 
previously  indicated  in  Fig.  9  (b),  were  used  to  weight  the  data  from  the  spectrometer.  This  relative  reflectance 
measurement  is  depicted  in  Fig.  10. 
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Fig.  10.  Relative  reflectance  of  the  healthy  and  unhealthy  leaves,  normalized  to  the  irradiance  in  order  2,  as 
measured  with  the  MSI  and  U2S.  The  healthy  vegetation  experiences  more  absorption  in  order  3  due  to  the 
presence  of  chlorophyll.  [11] 

The  root  mean  square  (RMS)  error  between  the  MSI  and  calibrated  spectrometer  data  is  calculated  by 


(15) 


The  RMS  error  for  the  healthy  and  unhealthy  leaf  samples  yields  0.0115  and  0.0117,  respectively.  Consequently,  the 
error  in  the  spectral  reconstruction  is  on  the  order  of  1 .2  percent  between  the  U2S  and  the  MSI.  A  2D  photo  of  the 
healthy  and  unhealthy  leaf  was  also  taken,  and  can  be  observed  in  Fig.  11.  The  healthy  leaf  has  a  lower  reflectance  in  the 
m  =  3  (red)  band  due  to  the  presence  of  chlorophyll. 
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Fig.  11.  Relative  reflectance  images  of  a  healthy  and  unhealthy  leaf,  (a)  Band-integrated  image  (m  =  0).  (b)  -  (e) 
Image  from  each  order.  m  =  2  through  m  =  5,  respectively.  [II] 

6.0  CONCLUSION 


The  development  of  a  snapshot  imaging  linear  polarimeter,  based  on  channeled  imaging  polarimetry  (C1P),  has  been 
overviewed.  By  incorporating  blazed  diffraction  gratings  into  each  arm  of  a  Sagnac  interferometer,  it  has  demonstrated 
that  white-light,  broadband  operation  of  the  C1P  technique  can  be  implemented.  Furthermore,  by  implementing  multiple- 
order  diffraction  gratings  into  the  system,  a  snapshot  multispectral  imager  can  be  realized.  Similar  to  how  C1P  modulates 
Stokes  parameters  onto  spatial  carrier  frequencies,  the  multispectral  implementation  of  the  sensor  produces  a  unique 
spatial  carrier  frequency  for  each  spectral  band. 
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